C语言解高次方程a*x^3+b*x^2+c*x+d=0

来源:百度知道 编辑:UC知道 时间:2024/06/08 14:43:04

刚才那个只能找一个根,现在改了下,可以找所有根:
刚才有个小BUG,改了下::
=================================================
#include <math.h>

void FindRoot(void)
{
double dbLeftX = -10.0;
double dbRightX = 10.0;
double dbResultY = 0.0;
double dbTmp = 0.0;
double dbMidleX =0.0;

//当二分点Y的值不为0不是根(即,二分点X不满足方程)
//当区间无限小时完成计算(下面的语句表示当区间不是无限小时继续计算)
while ( abs(dbRightX-dbLeftX) > 0.1e-6 )
{
//二分区间
dbMidleX = ( dbLeftX + dbRightX ) / 2.0;
//算二分点Y的值
dbResultY = 2 * pow(dbMidleX,3) - 4 * pow(dbMidleX,2) + 3 * dbMidleX - 6;

//如果计算的Y值为0时,找到一个根
if ( abs(dbResultY) < 0.1e-6 )
{
printf("Finde a root in [-10,10] , value = [%lf].\n",dbMidleX);
}

//判断方程的根落于刚才二分区间的哪个区间,并生成新的区间.
if ( dbResultY < 0.0 )
{
dbTmp = 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6;
if ( dbTmp < 0.0 )
{
dbLeftX = dbMid